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SUMMARY 

An analysis was made of the effect of centrifugal force on the 
flutter of a uniform cantilever beam. Two methods of analysis were 
used. The first method was an approximate method in which the 
assumption was made that the flutter motion could be represented by 
using Just two degrees of freedom. As these two degrees of freedom, 
the fundamental uncoupled modes in bending and torsion were used. 

In the second method, a solution .was obtained by solving the exact 
equations of motion of a beam fluttering in a centrifugal -force field, 
with no assumption made concerning the actual mode shapes at flutter. 
This solution is exact provided that the oscillatory aerodynamic 
forces are exactly known. A comparison of the results obtained by 
both methods showed a negligible difference, which indicated that 
the approximate method using two degrees of freedom could be used 
with safety. The results of the analysis Indicated that centrifugal 
force could be detrimental and could decrease the critical flutter 
speed. A beam that was stable at 0 rpm could become unstable at 
high rotative speeds. If, however, the flutter coefficient is kept 
below approximately 4, which is generally the case for compressor 
and turbine blades, there is little possibility of flutter occurring 
at low angles of attack. 


INTRODUCTION 

The study of vibrations of compressor and turbine blades has 
become of great importance as the use of axial-flow compressor and 
turbines has increased. Such studies involve resonant vibrations 
due to mechanical excitation from other parts of the machine or 
from pulsating air flows from previous rows of blades in addition 
to nonresonant aerodynamic excitation. The aerodynamic excitation 
involves a type of vibration that is self-sustained by the contin- 
ual absorption of energy from the air stream. This vibration is 
called flutter . 
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Classical flutter is more specifically defined as a self- 
sustained oscillation due to the coupling of inertia forces, elastic 
forces, damping forces, and dynamic aerodynamic forces. This type 
of flutter usually occurs on airplane wings at low angles of attack 
when the air velocity reaches a certain value called the critical 
flutter speed. It is to he distinguished from a radically different 
type of flutter called stalling flutter, which occurs on airfoils 
at high angles of attack, such as propeller blades near the stall 
point (reference l). Stalling flutter can be caused by an aero- 
dynamic hysteresis effect, the negative slope of the lift curve, or 
excitation by a system of Karman vortices (references 2 and 3). 

Classical-flutter theory has been developed by many investi- 
gators (references 4 to 6) and the results have been found to agree 
with experiment. A direct application of this unmodified classical 
theory to analysis of compressor and turbine blades would seem to 
indicate little possibility of such blades fluttering at zero 
angle of attack. This improbability of flutter is due to the stiff- 
ness of compressor and turbine blades compared with that of airplane 
wings. Two significant factors, however, must be considered: the 

effect of centrifugal force and the effect of cascading. These 
effects may be either beneficial, that is, increase the critical 
flutter speed or detrimental, that is, decrease the critical flutter 
speed . 


The classical approach to the flutter of a uniform wing is 
based on an important simplifying assumption, namely, the hypothesis 
of "semirigidity." This hypothesis implies that all flexural move- 
ments are in phase with one another and all torsional movements are 
in phase with one another. This hypothesis is equivalent to stating 
that the motion can be represented by a system of Just two degrees 
of freedom. For the two degrees of freedom, two arbitrary modes in 
bending and torsion are then chosen that are assumed to be inde- 
pendent of the air velocity. The fundamental uncoupled modes in 
bending and torsion of the wing vibrating in vacuo are usually 
assumed for these two modes. This assumption is made in order to 
make the problem mathematically tractable and has given good results. 

The exact solutions of the differential equations for the 
flutter of an elastic beam have been obtained (reference 7) for 
the special case of zero flexural stiffness and for the general 
case of finite flexural and torsional stiffnesses (reference 8). 

The results obtained show close agreement between the exact critical 
flutter speeds and the approximate values calculated on the basis 
of semirigid flutter theory. 
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An investigation was conducted at the NACA Lewis laboratory in 
order to obtain a rigorous solution to the problem of a uniform 
cantilever beam radially mounted fluttering in a cent if ugal -force 
field and to investigate the effect of centrifugal force on the 
flutter of such a beam. Most of the calculations are made by the 
approximate method using the assumption that the motion can be 
represented by the two fundamental uncoupled modes in bending and 
torsion. The fundamental torsional mode is assumed to be unaffected 
by the centrifugal force and the fundamental bending mode is cor- 
rected for the centrifugal force, as indicated by Timoshenko (refer- 
ence 9). The analysis then proceeds in the usual manner (refer- 
ence 4). The rigorous solution, which is exact if the oscillatory 
aerodynamic forces and moments are exactly known, is obtained by 
the use of matrix methods and is then used to determine the error 
involved in approximating the true flutter motion by the uncoupled 
modes in bending and torsion when centrifugal forces sure present. 


SYMBOLS 


The following symbols are used in the analysis : 


A(x) square matrix 

a coordinate of elastic axis measured from midchord, 

positive towards trailing edge, in units of half- 
chord 


a iJ 


elements of matrix 


a l' b l' c l, 

b 2 ,c 2 


functions of reduced frequency k and fundamental 
angular bending frequency 10 ^ or corrected funda- 
mental angular bending frequency co^'- 



l e> a e 


aerodynamic coefficients 


b y’ b y ,b e' b e 


c 

c 


half-chord used as reference unit length 
aerodynamic coefficients 

constant of integration 
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bending stiffness of beam 
nr 

function of k derived by Tbeodorsen (reference 4) 
functions of 

function of k derived by Theodorsen (reference 4) 

moment of Inertia about elastic axis per unit span 
length 

v=l 

torsional stiffness of beam 
reduced frequency, ~ 

oscillatory aerodynamic lift force per unit span 
length of beam 

oscillatory aerodynamic moment per unit span 
sum of moments acting on element of beam 
bending moment 
torsional moment 

mass of airfoil per unit span 

mass of cylinder of air (diameter of cylinder equal 
to airfoil chord) per unit span 

distance from center of rotation to root of beam 

location of center of gravity of airfoil measured from 
a, positive toward trailing edge in units of half- 
chord 

radius of gyration referred to a, in units of half- 
chord 


HACA TN 1693 


5 


S 

T 

t 

V 

y 

v/co^b 

x 


static mass unbalance of airfoil, (m)(r)(b) 

axial tension 

time 

shear force 

critical flutter speed 
flutter coefficient 

X * - P 

dimensionless distance along length of beam, — 


x' 


distance of element of beam from center of rotation 


Y(x) 

Y(0) 

y 

y 0 

y^x) 

y 2 (x) 

y 3’ y 4> y 5’ y 6 

6 


e 




column matrix of unknowns 

column matrix of constants 

vertical displacement 

maTlnuim bending amplitude 

dimensionless bending mode 

dimensionless torsional mode 

derivatives of y-^ and y 2 defined In text 

proportionality constant, function of radial distance 
from center of rotation to root of blade and of 
blade taper 

phase angle by which torsion vibrational mode lags 
bending mode during flutter 

angle of torsional displacement 

TnATimnm torsional amplitude 

ratio of mass of airfoil per unit span to mass of sur 

m 


rounding air cylinder per unit span. 
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square matrix denoting matrlzant of A(x) 
elements of 0 

frequency of flutter vibration 

fundamental bending frequency of airfoil 

fundamental bending frequency of airfoil corrected 
for centrifugal-force effect 

rotational frequency of beam 

fundamental torsional frequency of airfoil 


ANALYSIS 

The analysis of the effect of centrifugal force on the flutter 
of a uniform cantilever beam consists of two parts. The first part 
of the analysis uses the approximate solution to the flutter problem 
by the assumption of two degrees of freedom, namely, the uncoupled 
fundamental modes In bending and torsion. The fundamental bending 
mode is corrected for the centrifugal-force effect, whereas the 
fundamental torsional mode Is assumed to be negligibly affected by 
centrifugal force. Most of the results presented herein were 
obtained by this method. In the second part of the analysis, the 
validity of the approximate solution is determined by solving the 
exact equations of motion for a beam fluttering In a centrifugal- 
force field. The results are then compared with those obtained by 
the approximate method. In the entire analysis, the oscillatory 
air forces and moments acting an the beam are assumed to be the 
same as those acting on a two-dimensional airfoil In a uniform air 
stream. Incompressible flow is assumed and structural frictions 
are assumed negligible. 


Approximate Solution Assuming Two Degrees of Freedom 

Flutter equations. - The approximate equations for dynamic 
equilibrium of the system using the uncoupled fundamental modes In 
bending and torsion are given In reference 4. In the notation used 
herein they become 
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where 



where F and G are functions of only the reduced frequency k 
given in reference 4. 
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(4) 


then 


.. ,.c 

y = _ 0) y 


e = - co e 


(5) 


By combining equations (1) , (2), and (5) and dividing through 
by m a G £ 


M 


a y + la y + ' 11 


y + (v + la e - e ) 




be = o 


\ I 

1 

r»- 

co 

3 jby + 

|> - 1 »9 + - yj 


> (6) 


b 2 e = o 




Equations (6) have solutions other than y = 9 = 0 if and 
only if the determinant of the coefficients of y and 9 vanishes. 
If this determinant is set equal to zero, real and imaginary parts 
are separated, and the equation is rearranged. 


and 


M 3 

vt/ c 2 




a ll^t 


+ °1 ' 0 


J 


(7) 


b 2 = ^y + b e 


c 2 = ca y + U b 0 + a 0 (by - e) + b y (a 0 - e) - b^ - aj), 


J — Ok w _ 

ye 


\ - V <•* ' V + Ve - Vy ' (a 9 - a) (b y ' e) 


where 
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For a given airfoil, a 1 ,’b 1 ,c 1 ,‘b 2 ^ and Cg are functions only 

of k and the fundamental tending frequency corrected for centrif- 
ugal force effects 0)^', which is, in turn, a function of the rota- 
tive speed. For a given rotative speed, a, b, and c are func- 
tions of only the reduced frequency. Equations (7) can then be 
graphically solved for k and the ratio of flutter frequency to 
fundamental torsional frequency In the calculations per- 

formed herein, a value was assumed for k and the corresponding 
values of (o/co t and CO b 'A> t VQre calculated. This procedure 

obviated the necessity for a graphical solution. 

Bending frequency. - The effect of centrifugal force on fun- 
damental bending frequency 0)^ was obtained from reference 9. 

( O ^') 2 = ( i ) b 2 + P 2 ^ 2 ( Q ) 


where to b is the fundamental bending frequency without rotation 

p 1 b a function of the ratio of blade length to radial dis- 
tance from center of rotation to root of blade and of blade taper. 
By combining equations (7) and (8), k and the flutter fre^cy 
0) can be obtained as a function of the rotative speed. The crit- 
ical flutter speed v can then be obtained from the relation 


v = “5 (9) 

k 

Actually, it was thought more useful to obtain the dimensionless 
parameter v/oc^b as a function of the dimensionless parameter 

0> r /<V 
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Exact Solution for Beam Fluttering in 
Centrifugal-Force Field 

Equations of dynamic equilibrium. - The equations of motion 
for the beam are obtained from the equations for the equilibrium of 
forces acting on an infinitesimal element of the beam. 




The external forces and moments acting on an element of the 
beam as shown consist of transverse shearing forces and torsional 
moments that are transmitted from one portion of the beam to the 
next, plus the centrifugal force and moment and the aerodynamic 
force and moment. For dynamic equilibrium to exist, these external 
forces and moments must be balanced by the inertia forces and 
moments. The following equations of equilibrium can therefore be 
given: 
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For the equilibrium of moments in the x-y plane, 

M’ = ~ dx' - Vdx' - + T j-X- dx' + J2L L& l l = c 

dx' dx' 2 dx' ox' dx 1 2 

which upon neglecting terms of higher-order differentials becomes 

( 10 ) 


— - V + T <1^ = 0 
^irt ox' 


or 


Sx 

For the equilibrium of forces in the x' direction, 

2 

^7 dx' + m(o r x' dx' = 0 


- - m(0 r 2 i: 


which upon integration becomes 


mo) r 

2 

(*•) + 

T = ' 2 

when x' = 

Z + R, 

2 

m co r 

(Z + R) 

c = 2 


therefore 


T = 


m(i)J 2 2 


B) - (*') 


For the equilibrium of forces in the y direction. 


dV 


^7 dx' + L dx' - my dx' - S0 dx' =0 


( 11 ) 


( 12 ) 
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or 


dv 


^ t ' + L - my - S0 = 0 


From equation (10), 


. _ 2 + T " J + «S_ °Z_ 

^ atx'i 2 d(i ’) 2 




,2 

d y 


dr ch 


d(x')* &(*•)* 

If equations (ll), (12), (13), and (14) are combined, 


2 2 
o Mv m(0 _ 

£ + L_ 

d(x -) 2 2 


(Z + R ) 2 - (x ') 2 


n 2 

±JL 


d(x')‘ 


mo)- 2 x'^ + L- my-S0 = O 
r ox 


For the equilibrium of torsional moments, 


ST + M - SWj. 2 x- |L. - xe -Sy - 0 


From elementary elastic theory, 

d 2 ^ 


*<*') 


2 = - EI b 


d*7 

d(x -) 4 


and 


*t E d% 

^ d(x ') 2 


1893 


(13) 


(14) 


(15) 


(16) 


(17) 


equations (15) and (16) now become 


I 
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4 „ 2 

r 


cy* v m (0 


v * se + EI » suo 1 ' ~ 


(» + E ) 2 - ( x ’) 2 


A 


*2 

<Ll _ + 

Ja(x -) 2 


2 ,dy _ 
mt0 r X 3F " L = 


Sy + SO) 2 *•&- + ie - K- ■ ■— ^ - ? - M = 0 

r dx* d(x') 


(18) 




When the values of L and M as given by equation (2) are sub- 
stituted , 


my + S0 + EI^ — -j - 
d(x’) 


i- - (Z + R ) 2 - (x ') 2 




d(x') 


mCO 


_ n 2 2 

r + wm a (a y + la y) y + W m a b ( a fl +ia 0 ) e = 0 


a' y y 


0 0' 


2 2 

Sy + SO) ^x'l 2 - + 10 - K -- e ? + CO m b(b + ib )y + 

J r dx' d( x ') Z a y y 


•\ 


/ (19) 


2 2 — 

(0 m a b (b g + ib 0 )0 = 0 


J 


The following substitution of variables is convenient: 

x' - R 


x = 


( 20 ) 


Then equation (19) becomes 
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my 




8x' 


(1 + |) 2 - (x + |) 2 I 


dx 




mCO^ 2 (x-3) + « 2 m^a^ + i a v )y + 0) 2 m a b(a e +i a Q )0 = 0 


i ax a y 


and 


2 

Sy + SCO 2 (x4) |j + 10 - JL i-| + t/m b(b + 1 b )y 

r 1 ** i 2 bx 2 a 7 7 


CO 2 m a b 2 (b 0 + i b 0 )0 = 0 


>( 21 ) 


J 


These are the equations of motion. 

Solution of equations of motion. 

Let 

y = by 1 (x)e 


i(Ot 




and let 


then 


0 = yo'(x)e 


y 2 (x) = e" i8 y 2 ’(x) 


i( (Ot - 5) 


0 * y 2 (x)e 


ICOt 


J 


( 22 ) 


If these expressions are substituted Into equations (21) and 

2 loot 

the first equation divided by (0 m ft be and the second equation 

by O) 2 m a b 2 e i the following dimensionless equations are 

obtained : 
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(a y -n) + ia y 


EI b d4 yi liO) 2 


y + 1 —■ — — . - . ■■ JEL 


dx 20 ) 


p 2 p 2 

(l+j) -(x+|) 


,2 'N 

iii + 

djt 2 


, r. d yi 

2— (x + — ) + 

r\ ' 1 i 


0 )‘ 


l dx 


<a e - ur) + la e 


y 2 “° 


tv 


nr) + lb. 


(0 


B. d ^l 


7, + 4r — j- (* + 7 ) JJ- - 


K til 

< 0 2 m a ! 2 b 2 dx 2 


(O' 


= 0 


>(23) 


<V»*g ) + ib 0 


V 


J 


The first of equations (23) is a fourth-order linear differ- 
ential equation with variable coefficients and the second equation 
is a second-order linear equation with variable coefficients. No 
simple methods exist for solving such a pair of simultaneous equa- 
tions. The method used herein is to reduce the equations to a 
system of first-order equations and then apply the Peano-Baker 
method of integration (reference 10 ) . 


Let 




= £73 
dx “ y 4 


d 3 y x dy. 

^■ = 3T = y e 

dy 2 

dx” = y 5 




J 


(24) 


Then equations (23) became 
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where 


dy g 

dx” = a 61 y l + a 62 y 2 + a 63 y 3 + a 64 y 4 
dy 5 

d T = a 51 y l + a 52 y 2 + a 53 y 3 


~\ 


J 


(25) 


61 


4 - a - ia 

y y 

EIv 


03 2 m a l 4 


_ by - (ar + iby 

851 ” K 

0) 2 m aZ 2 h 2 


62 


nr - a e - Igg 
EI h 


(o^z 4 


, + lb i 

52 " K 


a_„ = 


- nd> r £ 


63 EI> 


(**!) 


uro) 


53 K 


M) 


m a Z 


_ 7 2v2 

V b 


H (0^ 
a 64 = 2EI b 

m Z 4 
a 


2 _ 


(1 + -|) - (1 + 7) 


B 2 


The system of six first-order equations (24) and (25) can now 
be replaced by the single first-order matrix equation: 


d - Y - M = A(x) Y(x) 
dx 


(26) 


where Y(x) is the column matrix 
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Y(x) = 


L y 6 


and A(x) is the square matrix: 



"0 

0 

1 

0 

0 0 


0 

0 

0 

0 

1 0 


0 

0 

0 

1 

0 0 

,A(x) = 

0 

O' 

0 

0 

0 1 


a 51 

a 52 

a 53 

0 

0 0 


_ a 61 

a 62 

a 63 

a 64 

0 0 

elements a. 

, of 

the 

matrix 

A(x) 

are not 


The solution of matrix equation (26) can now he set down as 


Y(x) = Y (0) 


(27) 


where is the matrizant of A(x) and is given hy 


J 
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n x = i + p x A dx + r aP x i a dx, dx + p x a r x i A p A d x 2 dx . d X 
o Jo Jo Jo 1 Jo Jo Jo 21 

. . . +P X A | 1 A. . . P* 11 A dx dx_ , . . . dx, dx + . . . 

Jo Jo Jo n 11-1 1 

00 

+ n ? 0 r C 1 a - • -r "' 1 C" a ai » 4vi • • • ax i ■ 


= i 


where x Q = x and I 1 b the unit matrix of order six. 

Boundary conditions. - The boundary conditions to be satisfied 
for the case of a cantilever beam are 


x = 0 (0) = y g (0) = y 3 (0) = 0 

x = 1 y 4 (1) = y g (1) = y g (1) = 0 


A 


(29) 


From the conditions at x = 0, the column matrix Y (0) becomes 

0 


Y (0) = 


0 

0 

J A (0) 

4(o. 

y 6 (o) 


For x = 1, equation (27) gives 

y 4 d) - 0 . Q 44 y 4 (0) ■+ 0 45 y 5 (0) + n 46 y g (o) 

y 5 (l) - 0 • n 54 y 4 <°> + Q 55 y 5 (0) + “56 y 6 (0) 

y 6 (1) - 0 - %i * 4 (°) * y.(o) + 0„„ 7.(0) 


(30) 




65 5 


66 6 


J 


(31) 
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where ft . . Is the element of row i and column j of the 
1J 1 

matrizant ft 


Critical flutter apeed and frequency, - The critical flutter 
speed and frequency can now he obtained by means of equations (31). 
Equations (31) have solutions other than y^(0) = yg(0) = yg(0) = 0 

if and only if the determinant of the coefficients vanishes. 
Therefore, 





°44 

°45 

n 46 








°S4 

n 5S 

°56 

= 0 



(32) 




°64 

^65 

°66 





or 










0 44 

n 55 

^55 

^64 

n 66 + 
- °4S 

^45 

°54 

^56 
^66 - 

^64 + 
°44 

^46 

^56 

^54 

^65 

^65 - 
= 0 

^46 

(33) 


The elements 



are complex and for a given airfoil at a 


given rotative speed are functions only of 0) and k. The real and 
imaginary parts of equation (33) can then be separated and the 
resultant two equations simultaneously solved for 00 and k. The 
calculations are very laborious and have been made herein for a 
given airfoil for two values of rotative speed. 


• RESULTS AND DISCUSSION 

Equations (7) were solved for k and (i)/co^ as functions 

of the ratio of torsional frequency to bending frequency corrected 
for centrifugal -force effect (O t /o) t 'j the results are shown in 

figure 1. Calculations were made for values of center- of- gravity 
position r of 0.1, 0.2, and 0.3. These values cover the prac- 
tical range. At r = 0.1, mass-ratio values p of 100, 300, 

700, and 1000 were used. At the two other values of r, values 
of p of 700 and 1000 were used. The calculations were made at 
such high values of p because the mass ratios of compressor and 
turbine blades are much higher than for airplane wings. The value 
of the radius of gyration r was taken as 0.5 for all calculations 
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and the value of elastic-axis position measured from midchord a 
as -0.4. Changes in a have a comparatively small effect (if the 
elastic axis is not very close to the gravity axis) and this one 
value in the practical range was chosen to reduce the number of 
calculations . 


Examination of figure 1 shows that the effect of incrasing 
U is to decrease both co/to ^ and k, the decrease in k 

generally being greater. Inasmuch as v is proportional to the 
ratio of oo / to to k, it will generally increase with increasing 

p. This effect becomes smaller at larger values of n. 

The effect of increasing r is to decrease the value of 
to / (0 ^ for high values of (o^/oo^ 1 , whereas the values 

of to /to t are left relatively unchanged at the low values of 

(0 /(O'. The values of k are either unaffected or increase at 

"t D 

the higher values of (0 1 /(0 ^ ' but are decreased at the lower 

values of (0 t /0)^' where the peak occurs. This relation would 

indicate that the effect of increasing r is to decrease v except 
near the point where (0^/u ' equals 1, where the peak occurs in 

the k curve . 


In order to determine the effect of centrifugal force on v, 
figure 1 must be more closely examined. Low values of the abscissa 
to^/to^' correspond to high values of rotative speed and high 

values of the abscissa correspond to low values of the rotative 
speed. This relation can be seen if equation (8) is rewritten in 
the form: 

a 

/ M tf ' 10 1 2 

V“ J ' *'• V (o 2 

d r 



1 + 3 



2 


(34) 
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For a given airfoil, 0) t /o)^' will therefore decrease as 

OO r /oo ^ increases. From figure 1 for values of co ^./co ^ ' Jess 

than approximately 1, the value of k will decrease and the value 
of Go/co^ will increase as the value of the rotative frequency 

0) increases (or as 0) /o> _ ' decreases). At values of 
r t b 

0) /CD • greater than approximately 1, both k and 00/(0 

D t 

increase as (0 increases (or as 00. /go' decreases). If k 
I* t D 

increases at a faster rate than 0) /go ^ in any part of the region 
in which both ( 0/(0 and k are increasing, a minimum flutter 

speed would exist at a certain rotative speed. This relation can 
be seen as follows: The critical flutter speed can be expressed as 



(35) 


Values of v will therefore increase with increasing rotational 
speed for values of less than approximately 1 and will decrease 
with increasing rotative speed for values of 00 ^/(o ^ ' greater 

than approximately 1. The existence of a minimum v is therefore 
indicated. 

The foregoing relations can perhaps be more clearly seen if 

equation (34) is used with figure 1 to plot directly k and 

(0 / CO as a function of the rotative speed. This relation is 

% 

shown plotted in figures 2 and 3 for values of <0 ranging 

from 1 to 10. The value of p is 700 and the constant 3 is 

set at 1.0, which corresponds to a nontapered blade with a ratio 

of root radius to blade length R/l of 0. As seen from figure 3 

as ( 0/(0 increases or as the rotative speed increases, the 
r' b 

frequency ratio go /go increases and from figure 2 the reduced 

"t 

frequency k in general at first increases and then decreases. 

If the rate of increase of k in any part of the region where it 
is increasing is greater than the rate of increase of (0 A°t' 

then from equation (35) a minimum value of v can be expected. 

In order to illustrate the previous conclusions, v was 
calculated from figure 1 and equation (35) for different values of 
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0) ^/o)^ and different values of R/l . The results are shown in 

dimensionless form in figures 4 to 6. The abscissa is (i) r /a) ^ 

and the ordinate is the dimensionless flutter coefficient v/w^b. 

The values of R/Z used are 0, 4 , and 10. The corresponding 
values of 6 taken from reference 9 are 1.00, 2.62, and 3.95. 

The values of v/ W^b, and for a given blade v, decrease 

to a minimum as co^/oo ^ or the rotative speed increases (figs. 4 

to 6). For a given value of u and r, this minimum value of 
v/w^b is a constant for all values of OJ^/co^. The value of 

W r /(0^ at which this minimum occurs, however, increases as the 

ratio CO ^/co ^ increases . This decrease in v as the rotative 

speed is increased is caused by the centrifugal force causing only 

CO , to increase . Although increasing both CO and (0 A will 
b b t 

increase v, increasing CO^ alone may decrease v (figs. 4 to 

6 and reference 11). Values of v may be decreased by the 
centrifugal-force effect as much as 65 percent. 

The effect of R /Z is in the same direction as the effect of 
rotative speed. Increasing R/l increases the centrifugal-force 
effect. The effect is therefore the same as increasing the rotative 
speed. The minimum points on the curves are therefore shifted to 
lower values of the rotative speed as can be seen in the figures. 

The lowest value of the flutter coefficient v/co^b obtained 

in these calculations was approximately 4.3. In the practical 
range of parameters for current compressor and turbine blades, this 
value is probably close to the lowest value that is to be expected. 
If the blades are so designed that the operating v/ to t b is lower 

than approximately 4, little possibility exists of classical flutter 
occurring. When sufficient knowledge becomes available on the 
effects of cascading on the oscillatory aerodynamic forces and 
moments, the previous results may have to be modified. The previous 
results apply to low angles of attack. With blades operating at 
high angles of attack near the stall point, v may be considerably 
reduced . 

In order to determine the possible error in the previous cal- 
culations because of approximating the flutter motion by the assump- 
tion of uncoupled modes in bending and torsion, an exact solution 
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was obtained for two points by means of equation (33) . The results 
are shown in figure 7 where k and w/io^ are plotted against 

to^/to-^. The curves were obtained by the approximate method used 

in obtaining figures 1 to 6. The circles represent the points 
obtained by the exact solution. The difference in the results 
obtained by the approximate solution and by the much more laborious 
exact method is negligible . The point at CO / to ^ = 0 verifies 

the conclusions of reference 8 that the use of the uncoupled modes 

in bending and torsion as degrees of freedom gives results that 

differ negligibly from the exact solution for a simple cantilever 

beam. The point at to /to. = 1 indicates that the same conclu- 

r' b 

sion can be drawn for a cantilever beam in a centrifugal- force field 
provided that the fundamental bending mode is corrected for 
centrifugal-force effect . 


CONCLUSIONS 

On the basis of calculations carried out by exact and approxi- 
mate methods on the effect of centrifugal force on the flutter of a 
uniform cantilever beam it is concluded that : 

1. The approximate method of solution, using as two degrees 
of freedom the uncoupled fundamental modes in bending and torsion 
in vacuo, gave results that differed negligibly from those obtained 
by exact methods, which indicated that the approximate method can 
be used with safety. 

2. The effect of centrifugal force could be detrimental and 
could decrease the critical flutter speed for the range of parameters 
considered as much as 65 percent under certain conditions. A beam 
that was stable at 0 rpm could become unstable at high rotative 
speeds. 

3. In the practical range of blade parameters, ignoring the 
effect of cascading, there would be little probability of compressor 
and turbine blades fluttering at low angles of attack if the flutter 
coefficient was kept below approximately 4. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 6, 1949. 
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(a) Center-of -gravity position measured from position 
a r, 0,1; elastic-axis position measured from 
midchord a, -0.4; radius of gyration r g , 0.5. 

Figure 1. - Variation of reduced frequency and ratio of 
flutter frequency to fundamental torsional frequency 
with ratio of fundamental torsional frequency to 
corrected fundamental bending frequency for various 
mass ratios. 
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(b) Center-of -gravity position measured from position 
•a, r, 0.2; elastic-axis position measured from 
midchord a, -0.4; radius of gyration r g , 0.5. 

Figure 1. .- Continued. Variation of reduced frequency 
and ratio of flutter frequency to fundamental torsional 
frequency with ratio of fundamental torsional frequency 
to corrected fundamental bending frequency for various 
mass ratios. 
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(c) Center-of-gravity position measured from position 
a, r, 0,3; elastic-axis position measured from 
midchord a, -0.4; radius of gyration r g , 0.5. 

Figure 1. - Concluded. Variation of reduced ' frequency 
and ratio of flutter frequency to fundamental torsional 
frequency with ratio of fundamental torsional frequency 
to corrected fundamental bending frequency for various 
mass ratios. 
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Rotational frequency t Qr 

Fundamental bending frequency 

Figure 2. - Variation of reduced frequency with ratio of 
rotational frequency to fundamental bending frequency 
for various ratios of fundamental torsional frequency 
to fundamental bending frequency. Mass ratio 700; 
center-of -gravity position measured from position a, r, 
0.2; root-radius-to-length ratio R/l , 0. 
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Figure 3. - Variation of ratio of flutter frequency to 
fundamental torsional frequency with ratio of rotational 
frequency to fundamental bending frequency for various 
ratios of fundamental torsional frequency to fundamental 
bending frequency. Mass ratio n, 700; center-of -gravity 
position measured from position a, r, 0.2; root-radius- 
to-length ratio R/l, 0. 
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(a) Maas ratio 4* 300; root -radius -to -length 
ratio R/l 9 0. 



(b) Mass ratio 300; root-radius- to- length 
ratio R/l 9 4. 


Figure 4. - Variation of flutter coefficient with ratio 
of rotational frequency to fundamental bending frequency 
for various ratios of fundamental torsional frequency 
to fundamental bending frequency. Center-of -gravity 
position measured from position a, r,. 0.1; elastic- 
axis position measured from midchord a, -0.4; radius 
of gyration r , 0.5. 
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(d) Mass ratio H , 700; root-radius- to- length 
ratio R/l, 0* 

Figure 4. - Continued. Variation of flutter coefficient 
with ratio of rotational frequency to fundamental 
bending frequency for various ratios of fundamental 
torsional frequency to fundamental bending frequency. 
Center-of -gravity position measured from position a, r, 
0.1; elastic-axis position measured from midchord a, 
-0.4; radius of gyration r . 0.5. 

ft 
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(f) Mass ratio \i , 700; root -radius -to-length 
ratio R/l » 10. 


Figure 4. - Concluded. Variation of flutter coefficient 
with ratio of rotational frequency to fundamental 
bending frequency for various ratios of fundamental 
torsional frequency to fundamental bending frequency. 
Center-of-gravity position measured from position a, r, 
0.1; elastic-axis position measured from midchord a, 
-0.4; radius of gyration r , 0.5. 




Flutter coefficient 


NACA TN 1893 


33 




(b) Maas ratio M*, 700; root -radius -to -length 
ratio R/l , 4. 


Figure 5. - Variation of flutter coefficient with ratio 
of rotational frequency to fundamental bending frequency 
for various ratios of fundamental torsional frequency 
to fundamental bending frequency. Center-of -gravity 
position measured from position a, r, 0.2; elastic- 
axis position measured from midchord a, -0.4; radius 
of gyration r , 0.5. 
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Figure 5. - Concluded. Variation of flutter coefficient 
with ratio of rotational frequency to fundamental 
bending frequency for various ratios of fundamental 
torsional frequency to fundamental bending frequency. 
Center-of -gravity position measured from position a, r, 
0.2; elastic-axis position measured from midchord a, 
-0.4; radius of gyration r g, 0.5. 
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(a) Mass ratio M- , 700; root radius to length 
radius R A » 0. 
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Rotational frequency ^r 

Fundamental bending frequency' 

(b) Mass ratio 700; root-radius -to-length 
ratio R/l , 4. 

Figure 6. - Variation of flutter coefficient with ratio 
of rotational frequency to fundamental bending frequency 
for various ratios of fundamental torsional frequency 
to fundamental bending frequency. Center-of -gravity 
position measured from position a, r, 0.3; elastic- 
axis position measured from midchord a, -0.4; radius 
of gyration r , 0.5. 
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(c) Mass ratio ijl, 700; root-radius-to-length 
ratio R/l, 10. 

Figure 6. - Concluded. Variation of flutter coefficient 
with ratio of rotational frequency to fundamental 
bending frequency for various ratios of fundamental 
torsional frequency to fundamental bending frequency. 
Center-of-gravity position measured from position. a, r, 
0.3; elastic-axis position measured from midchord a, 
-0.4; radius of gyration r g , 0.5. 
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